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ABSTRACT 

We analyze the temporal variation of the diurnal anisotropy of sub-TeV cos- 
mic ray intensity observed with the Matsushiro (Japan) underground muon de- 
tector over two full solar activity cycles in 1985-2008. We find an anisotropy 
component in the solar diurnal anisotropy superimposed on the Compton-Getting 
anisotropy due to the earth's orbital motion around the sun. The phase of this 
additional anisotropy is almost constant at ~15:00 local solar time corresponding 
to the direction perpendicular to the average interplanetary magnetic field at the 
earth's orbit, while the amplitude varies between a maximum (0.043±0.002 %) 
and minimum (~0.008±0.002 %) in a clear correlation with the solar activity. We 
find a significant time lag between the temporal variations of the amplitude and 
the sunspot number and obtain the best correlation coefficient of +0.74 with the 
sunspot number delayed for 26 months. We suggest that this anisotropy might be 
interpreted in terms of the energy change due to the solar-wind induced electric 
field expected for GCRs crossing the wavy neutral sheet. The average amplitude 
of the sidereal diurnal variation over the entire period is 0.034±0.003 %, which 
is roughly one third of the amplitude reported from AS and deep-underground 
muon experiments monitoring multi-TeV GCR intensity suggesting a significant 
attenuation of the anisotropy due to the solar modulation. We find, on the other 
hand, only a weak positive correlation between the sidereal diurnal anisotropy 
and the solar activity cycle, in which the amplitude in the "active" solar ac- 
tivity epoch is about twice the amplitude in the "quiet" solar activity epoch. 
This implies that only one fourth of the total attenuation varies in correlation 
with the solar activity cycle and/or the solar magnetic cycle. We finally examine 
the temporal variation of the "single-band valley depth" (SBVD) quoted by the 
Milagro experiment and, by contrast with recent Milagro's report, we find no 
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steady increase in the Matsushiro observations in a 7-year period between 2000 
and 2007. We suggest, therefore, that the steady increase of the SBVD reported 
by the Milagro experiment is not caused by the decreasing solar modulation in 
the declining phase of the 23rd solar activity cycle. 

Subject headings: diurnal anisotropy of sub-TeV cosmic ray intensity — solar 
modulation of the sidereal anisotropy — cosmic ray observation with underground 
muon detector 
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Introduction 



The sidereal anisotropy of high-energy galactic cosmic ray (GCR) intensity provides 
us with unique information about the magnetic structure of the heliosphere and/or the 
local interstellar space surrounding the heliosphere, through which GCRs propagate to the 
earth. Examining the temporal variation of the anisotropy is also of interest, since such 
a variation, if any, may reflect the change in the global magnetic structure in response to 
the solar activity- and/or magnetic-cycle variations. While recent observations with air 
shower arrays and underground muon detectors have reported a consistent global structure 



of the anisotropy, cone 



N agashi ma et al. 



2009; 



1989 



usive r eports on its temporal variation are stil 



Amenomori et al. 



2004: 



Amenomori et al. 



2005 



Guillian et al. 



limited (e.g., see 

~3 



200 



Abdo et al. 



20091 ). The galactic anisotropy is measured as a sidereal daily 
variation (SDV) of cosmic-ray intensity recorded in a fixed directional channel of the 
detector on the spinning earth. The temporal variation of the SDV in connection with 
the solar magnetic cycle was investigated from a continuous observation of the SDV of 
~10 TeV cosmic-ray intensity over 12 years between 1973 and 1987 with an air-shower 
(AS) detector at Mt. Norikura in Japan. While the SDV in each year was statistically 
significant, i ts temporal variation wa s quite small with a statistical significance of only 



two sigmas ( jNagashima et al. 



19891 ) . Comparing data divided into two 5-year periods, 
1997-2001 including the solar maximum and 2001-2005 approaching the solar minimum of 
the 23rd solar activity cycle, the Tibet III air shower experiment also concluded that the 
variation of 3 TeV cosmi c ray anisotropy during the solar activity cycle was insignificant 



( lAmenomori et al. 



20061 ). In the present paper, we analyze the long-term variation of the 
SDV by analyzing 0.6 TeV cosmic ray data observed with a Japanese underground muon 
detector over 24 years between 1985 and 2008. This observation period is long enough to 



examine the variation of both the 11-year solar activity- and 22-year solar magnetic-cycles. 
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Recently, the Milagro experiment utilizing a large water Cherenkov detector 
surrounded by an air shower array, on the other hand, has reported on a steady increase in 



the m agnitude of t 



re an isotropy of 6 TeV cosmic ray intensity over 7 years between 2000 and 



2007 (AbdoetaL 



20091 ) . This period corresponds to the declining phase of the 23rd solar 
activity cycle when the yearly mean sunspot number decreased from ~100 to ~10 toward a 
long-lasting activity minimum epoch. The increase reported by the Milagro experiment was 
so large that the anisotropy amplitude more than doubled between 2000 and 2007. If this 
increase did indeed, result from the decreasing solar activity, it is reasonable to expect that 
the effect should be observable also in the sub-TeV region, because the physical processes 
responsible for the solar modulation, such as adiabatic deceleration, pitch angle scattering 
by the magnetic irregularities and the solar wind convection effects in the heliosphere, 
are all expected to be more effective for lower energy cosmic rays. It has actually been 
established that the amplitude of the sidereal anisotropy decreases with decreasing energy 
between ~0.1 and ~1 TeV due to a larger attenuation by the solar modulation effect, while 
the amplitude is more constant for multi-TeV cosmic rays with energie s over ~1 TeV, 



implying insignificant influence of the solar modu 



1997 



Nagashima et al. 



1998 



Amenomori et al. 



ation on Te V GCRs (IMunakata et al. 



2005 



2009al ). In the present paper, we 



also examine whether the effect reported by Milagro is observed in the sub- TeV cosmic ray 
anisotropy. 

The outline of this paper is as follows. In §2, we briefly describe the data and analysis 
method and present analysis results. We first show in §2.1 how the diurnal anisotropy in 
solar time changes in a significant correlation with the solar activity, indicating that solar 
modulation effects act on sub- TeV GCRs. We then show in §2.2 the long-term variation in 
the sidereal anisotropy. Our results are discussed in §3. 
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2. Data and analysis 

The Matsushiro underground muon detector (36.53°N, 138.01°E, ~300 km north-west 
of Tokyo) has been in operation since 1984. With its vertical overburden of 220 meters 
water equivalent, its muon energy threshold is 100 GeV. The multi-directional muon 
detector (sometimes also referred to as Zohzan according to the name of the hill above the 
observatory) consists of two horizontal layers of plastic scintillators, vertically separated 
by 1.5 m. Each layer comprises a 5x5 square array of 1 m 2 unit detectors, each with a 
1 mxl mxO.l m plastic scintillator viewed by two photomultiplier tubes (PMTs) of 12.7 
cm diameter. In order to reduce the influence of the dark current noises, each detector is 
designed to output signals only when two PMTs simultaneously yield pulses each shaped 
with a width of 300 nanosecond. By counting pulses of the twofold coincidences between 
a pair of detectors on the upper and lower layers, we record the rate of muons from the 
corresponding incident direction. The multidirectional muon telescope comprises various 
combinations between the upper and lower detectors. In this paper, we analyze hourly 
muon rates in the eastern (E) and western (W) viewing channels comprised of 17 directional 
channels available from the detector. The average count rates of E and W channels are 
7.14 Hz and 9.89 Hz, respectively, while the median primary cosmic ray energies recorded 
in these two composite channels are respectively calculated to be 0.660 TeV and 0.549 TeV 
on t he basis of the response function of the atmospheric muons to the primary cosmic 



ray (Murakami et al. 



19791 ). The median energy of primary cosmic rays recorded in the 



vertical channels is also calculated to be 0.596 TeV. 

In order to properly eliminate the spurious variations due to the atmospheric effect, 
we adopt the following "East- We st" method which is similar to the "Forward-Backward" 



method applied to Milagro data (lAbdo et al. 



20091 ). Let I(t) be the fractional deviation of 
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the intensity in a vertical directional channel at the time t in hours, defined as, 



I(t) = [N(t)-N)/N, (1) 

where N(t) is the hourly muon count and N is the 24-hour central moving average of N(t) 
at t. Also let R(t) be the percent deviation due to the anisotropy expected in the vertical 
directional channel. By taking the difference between I(t)s in the East- and West-viewing 
channels, I E (t) and I (t), we can deduce the "differential" variation of R(t) as, 



D(t) = dR(t)/dt = [I h (t) - I w (t)]/At 



(2) 



where At is the hour angle separation between the mean East- and West-incident directions 
averaged over the East- and West-inci dent muons and J E (t) and I w (t) are respectively 



composed of five directional channels (jMunakata et al. 



20061 ). By utilizing the response 



function of the atmospheric muons to p rimary cosmic ra ys together with the particle 



trajectory code in the geomagnetic field (jLin et al 



19951 ). we calculate At to be 5.3 hours. 



Note that the "differential" variation D(t) in eq. (2) is free from the spurious variation 
due to the atmospheric effects which is expected to be common for I E (t) and I w (t). We 
consider an expansion of R(t) into a Fourier series expressed with the amplitude and the 
phase 4>ni as 

3 

R{t) = ^a^cosnw(t-0^), 

n=l 

where u — ir/12. By introducing this R(t) into eq. (2), we obtain D(t), as 

3 



(3) 



D(t) = a E cos nu(t — 



(4) 



n=l 



where the coefficients a® and 4>^ are respectively related to and 0^, as 



D 



+ 6.0. 



(5) 
(6) 
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We deduce a® and every year by best-fitting eq. (4) to the observed yearly mean D(t), 



as 



tf = tajr\B2/A2), 



where 



24 



A® = (w/tt) d (U) cosuU, 
i=i 

24 

: (ou/ir) -D(tj) sin uti. 



B 



i=i 



Errors of and 0^ are propagated from the count rate error <5j of D(ti), as 



% = (cos^/a^)V(^^) 2 + (^^) 2 , 



where 



<L4 



(w/tt) 



24 



a y^ 2 cos 2 a;^, 
\ i=i 



^ = ( W /7T), 



24 



a ^ 2 sin 2 a;tj. 
^ i=i 



(7) 
(8) 



(9) 
(10) 



(11) 
(12) 



(13) 



(14) 



We then obtain a^, 0^ and errors according to eqs. (5) and (6). We also calculate the daily 
variation, i.e. the 24- hour profile of R(t), by "integrating" D(t) in eq. (2) with respect to t 
as 

i 

R*(t t ) = J2D(t k )-(R*), (15) 



k=i 



where (R*) = £?=i £!=i D{t k )/2A. 



Fig. 1 displays af L and 0^ for n — 1 as harmonic vectors in the sidereal (SI) (366.25 
days yr" 1 ), solar (SO) (365.25 days yr _1 ) and anti-sidereal (ASI) (364.25 days yr" 1 ) time 
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frames. The vectors in ASI time frame are nonphysical and are calculated for checking the 
"sideband" effect possibly arising from the seasonal variation in the amplitude of vectors in 
the SO time frame. It is clear that the diurnal (n = 1) vectors are most significant in SI 
and SO time frames, while they are less significant in ASI time frame ensuring the small 
sideband effect. The amplitude of the semi-diurnal (n = 2) vector is about one third of the 
diurnal vector and is significant in SI and SO time frames, while it is insignificant in the 
ASI time frame. The tri-diurnal (n = 3) vector is insignificant in all time frames. 

The observed diurnal vector in the ASI time frame in Fig. 1 is consistent with the 
sideband effect expected from the seasonal change of the diurnal v ector arising from 



the second order anisotropy of solar origin ( jNagashima et al. 



19831 ). The second order 



anisotropy is symmetric with respect to the interplanetary magnetic field with minimum 
intensities in directions parallel and anti-parallel to the magnetic field. This anisotropy 
produces the diurnal variation as well as the semi-diurnal variation in the SO time frame 
because of the inclination of the earth's spin axis from the magnetic field direction. As the 
angle between the earth's spin axis and the mean magnetic field is subject to the seasonal 
variation due to earth's orbital motion around the sun, the amplitude of the diurnal 
variation also shows the seasonal variation producing spurious diurnal variations in the ASI 
and SI time frames as sideband effect. In this case, the space harmon ic vectors in the AS I 



and SO time frames are expected to satisfy the following conditions (jNagashima et al. 



19831 ). 



^RiAsi) = 23 3 hour) (16) 

a n(Asi )/a R(so) = 3g) (lg) 

where c^f^ 57 ^ (4>2^ S °^) an d af (a^' SO ' ) ) denote respectively the phase and amplitude 
of the diurnal (semi-diurnal) vector in the ASI (SO) time frame in space. The observed 



- 10 - 



average (p^ ASI ^ {(p^ ^) i n the ASI (SO) time frame over the entire observation period 
is 2. Oil. 8 hour (2.5±0.6 hour), while the observed amplitude ratio af^ ASI ^ /a^ ^ is 
0.26±0.15, after the correction for the geomagnetic deflection and the attenuation in the 
atmosphere. The observed harmonic vectors, therefore, are consistent with the expectation 
within er rors. Based on thi s , we o btain the spurious harmonic vectors in SI time frame 



following iNagashima et al. 



(I1983 ). as 

a R(si) = .948af (AS7) , (19) 

^RiSI) = ^R(ASI) _ 4 53 hou ^ (2Q) 

a? 5/) = 0.098a? 5O) , (21) 
^(si) = ^(so) + 4 56 hQur _ (22) 

We correct the observed harmonic vectors in the SI time frame by subtracting spurious 
vectors in eqs. (19)-(22). Fig.l displays the corrected and uncorrected SI vectors by full 
black and full gray circles, respectively. Table 1 presents the corrected and 0f in the SI 
time frame, together with the observed af and 0f in the SO and ASI time frames for each 
year between 1985 and 2008. We analyze these results in the following subsections. 



2.1. Solar cycle variation of the solar diurnal anisotropy 



One remarkable feature of Fig. 1 is that the solar diurnal vector changes periodically 
every ~11 years. This is evident from Fig. 1 as the summation dial of the diurnal vector 
in the SO time frame is winding periodically in the second quadrant between 06:00 
and 12:00 local time (LT). In our separate paper, we reported that this is due to an 
"additional" anisotropy component superposed on the constant Co mpton- Getting (CG) 



anisotropy arising from t 



1935 



Amenomori et ai- 



re earth's orbital motion aro und the sun (jCompton fc Getting 



2004 



Munakata et al. 



20061 ). This "additional" anisotropy has 
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a constant phase at ~15:00 local time (LT) and an amplitude changing periodically every 
~11 years. The year-to-year variation of the phase and amplitude of this additional vector, 
derived by subtracting the expected CG effect from the observed vector, is shown in Fig. 2b. 
For comparison, we also show in Fig. 2a the yearly mean sunspot number (SSN) and the 
neutral sheet tilt angle (TA), which are often used as indicators of the solar activity and 
the solar modulation. The average phase of the additional anisotropy in the 6 years (1991, 
1992,1994 and 2002-2004; period I) when the amplitude exceeds 0.035 % is 15:10±00:25 
LT, while the average amplitude in the same period is 0.043±0.002 %, with errors deduced 
from the dispersion of 6 yearly values. On the other hand, the average amplitude in the 6 
years (1986-1988, 1997, 1998 and 2007; period II) when the amplitude is smaller than 0.015 
% is 0.008±0.002 %. Fig. 3 shows the average solar daily variations observed in period I 
(a) and period II (b). Plotted in this figure is R*(ti), obtained by "integrating" D(t) in 
eq. (15). It is clear that R*(ti) is fairly consistent with the CG anisotropy in panel (b) with 
the maximum intensity at ~06:00 LT displayed by a dashed curve, while R*(ti) in panel (a) 
deviates from the CG anisotropy showing a maximum intensity at ~12:00 LT due to the 
additional anisotropy with the maximum at ~15:00 LT superimposed. 

While a clear ~ 11-year solar cycle variation is seen in the amplitude of the additional 
anisotropy in Fig. 2b, there is a significant time lag between the temporal variations of 
the amplitude and the SSN or TA in Fig. 2a. We calculate the 12-month central moving 
average of the SSN and evaluate the correlation coefficient between the amplitude and the 
SSN shifting the SSN data to the later period every one month from to +48 months. 
The time-shifted coefficient gradually increases from +0.25 at month lag, reaching the 
maximum of +0.74 at +26 month lag and then decreases to +0.26 at +48 month lag. We 
also calculate the correlation with the TA plotted in Fig. 2a and obtain the maximum 
correlation coefficient of +0.63 at +30 months. Based on these observational results, we 
will discuss the possible origin of the additional anisotropy in §3. 
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2.2. Long-term variation of the sidereal anisotropy 



Fig. 2c shows the year-to-year variation of the amplitude and phase of the sidereal 
diurnal anisotropy. The average amplitude over the entire period is 0.034±0.003 % 



muon experiments monitoring multi-TeV GC : 


I intensity (e.g., see 


Naeashima et al. 


1989: 


Amenomori et al. 


2005; 


Guillian et al. 


2007; 


Abdo et al. 


2009 




Amenomori et al. 


2009) 



This attenuation is probably due to the solar modulation effects, but there is only a weak 
positive correlation between the amplitude and the SSN and/or TA in this figure. Such a 
weak positive correlation also has been reported f rom the long-term obser vation using an 



underground muon detector at a shallower depth ( jNagashima et al. 1 120041 ). We obtain the 
maximum correlation coefficient of +0.57 (+0.50) between the amplitude and the SSN (TA) 
at the time lag of +14 (+14) months applied to the SSN (TA). We examine the correlation 
with the solar activity by comparing the average sidereal daily variations R*(ti) in eq. (15) 
in the "active" period of 9 years (1988-1992, 1999-2002) when the yearly mean SSN exceeded 
90 with that in the "quiet" period of 11 years (1985-1987, 1994-1997, 2005-2008) when the 
SSN was below 30. In order to examine the correlation with the 22-year solar magnetic cycle, 
we also calculate the average R*(ti) in the "positive" (8 years; 1991-1998) and "negative" 
(12 years; 1985-1988, 2001-2008) periods, separated by the polarity reversal epochs of the 
solar polar magnetic field indicated in Fig. 2. During the "positive" ("negative") period, 
the average interplanetary magnetic field (IMF) is directed away from (toward) the sun in 
the northern hemisphere. Fig. 4a-4d display R*(ti)s in the "active", "quiet", "positive" 
and "negative" periods, respectively. The differences between R*(ti)s in the four periods 
in Fig. 4 are all insignificant. The amplitude of the daily variation in "active" period in 
Fig. 4a, for instance, appears larger than that in the "quiet" period in Fig. 4b, but the 
intensity difference in each LT remains below three times of the error deduced from the 
dispersion of yearly mean intensities in each period. The difference between R*(ti)s in 



13 



"positive" and "negative" periods in Figs. 4c and 4d remains below two times of the error. 
It is evident, on the other hand, that the overall feature of R*(ti) in Fig. 4 with a minimum 
intensity at ~12:00 LT and a maximum at ~06:00 LT is common for all periods and is 
quite similar to that obser ved with AS and deep-underground muon experiments measuring 



multi-TeV GCRs (e.g., s ee 



2007 



Abdo et al. 112009 



Nagash ima et al 



Amenomori et al. 



1989 



20091 1 . 



Amenomori et al. 



2005 



Guillian et al. 



We finally compare the temporal variation of R(t) with that reported by the Milagro 



experiment. Fig. 5 c ompares the 



Milagro experiment (I Abdo et al. 



'sing 



e-band valley depths" (SBVDs) reported by the 
20091 ) and those observed by Matsushiro during the 
identical 7-year period from 2000 to 2007. The SBVD by Matsushiro in this figure is 
multiplied by three in order to roughly compensate the attenuation of the amplitude in 
sub-TeV region. We derived the SBVD from Matsushiro data by calculating the minimum 



intensity of R(t) reconstructec 



analysis period in 



Abdo et al. 



rom yearly mean a® and (p® i n eq. (3). Adjusting to the 



(2009), we calculated the yearly mean SBVD by Matsushiro 



from 12-month data between June and July. The steady increase in the SBVD reported 
by the Milagro experiment is not seen in the Matsushiro record. We also calculated SBVD 
every year in the entire period of the observation by Matsushiro and confirmed that the 
long-term variation of the SBVD is similar to the variation of the amplitude of the sidereal 
diurnal anisotropy in Fig. 2c which shows no significant correlation with the solar activity- 
and magnetic-cycles. We conclude, therefore, that the steady increase of the SBVD reported 
by Milagro experiment is, most likely, not due to the decreasing solar modulation in the 
declining phase of the 23rd solar activity cycle. 
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3. Discussions 



A clear ~1 1-year variation was found in the solar diurnal anisotropy. The amplitude 
of the residual anisotropy, after subtracting the contribution from the Compton-Getting 
anisotropy due to the earth's orbital motion around the sun, varies in an ~ll-year cycle 
between a maximum amplitude of 0.043±0.002 % and a minimum of 0.008±0.002 %, 
while its phase remains fairly constant around 15:10±00:25 LT. A possibility is that the 



anisotropy is connected with the temporal vari ation of the IM 



Erdos&Kota 



whic h, on large scales, is 



(Il980l ) first suggested that 



influenced by the neutral sheet tilt angle (TA). 
GCRs are expected to experience the energy change when they cross the neutral sheet 
on their orbital motion through the heliosphere to the earth. They computed the energy 
change along individual trajectories and deduced the directional anisotropy expected at the 
earth. Preliminary results of recent numerical model simulations using different tilt angles 



( Kota et al. 



20071 ) showed some qualitative similarities to the Matsushiro observations for 
solar anisotropies: high tilt angles in the maximum period produced larger anisotropies 
while the phases turned out fairly stable typically in the 3rd quadrant between 12:00 and 
18:00 LT. We tentatively note that the ~15:00 phase corresponds to arrival direction of 
particles that have the best chance to gain energy from the interacting with the corotating 
magnetic sectors divided by a wavy current sheet. It is most interesting, on the other hand, 
that the variation of the observed amplitude of the "additional" anisotropy, exhibits a 
significant time lag of +26 months relative to the SSN. This time lag corresponds to the 
time for the solar wind with an average speed of 400 km/s to propagate through ~180 AU, 
which is larger tha n the heliocentric distance of the solar wind termination shock (TS) in 



the nose direction (IDecker et al. 



2005 



Stone et al. 



20081 ). Even with the slowing down of 



the subsonic solar wind, this ti me-lag corresponds to dista nces deep in the heliosheath, and 



possibly close to the heliopause (IWashimi fc Tanaka 



1996J). The cyclic variation induced by 



the sun propagates to the outer heliosphere with the average solar wind speed. The global 
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distribution of the solar wind plasma and the magnetic field, which governs the TA and the 
spatial distribution of GCRs, is reorganized when the variation propagates throughout the 
entire region filled with the solar wind. This implies, therefore, that the intensity of 0.6 TeV 
GCRs is still a subject to the solar modulation operating over an entire region within the 
TS and possibly also in the subsonic region beyond the TS. The possible role of heliosheath 
outside the TS is largely unknown at present and remains to be explored. 

We discuss the sidereal anisotropy, which shows no clear correlation with either the 
solar activity- or magnetic-cycles. Furthermore, there is even a factor of three attenuation 
of the amplitude relative to air shower observations due to the solar modulation. A possible 
explanation for this may be found in the region where the attenuation takes place. The solar 
diurnal anisotropy of 0.6 TeV GCRs at the earth's orbit is mainly produced by the global 
distribution of the GCR density inside the TS, while the sidereal anisotropy originates from 
the GCR anisotropy in interstellar space outside the heliosphere and it is attenuated as 
cosmic rays propagate to the earth through the heliosphere losing their original directional 
information. If the attenuation occurs predominantly inside the TS, it is reasonable to 
expect significant solar cycle variation of the anisotropy at the earth's orbit, as we saw in 
the solar diurnal anisotropy. The situation may change, however, if the attenuation occurs 
in the region outside the TS. Although our knowledge about the heliosheath between the 
TS and the heliopause is still limited, the Magneto-Hydrodynamic (MHD) simulations of 
the interaction between the solar wind and the interstellar plasma suggest the existence of a 



lon g heliotail extending over thousand of AU downstream the interste 



see 


Baranov & Malama 


1993 


1998: 


Pogorelov et al. 


2004 


)■ 



Washimi fc Tanaka 



1996 



Zank et al. 



ar p 



1996 



asma flow (e.g ., 



Linde et al. 



the larger extent of the heliotail, the propagation time of the solar cycle variation through 
the entire heliotail is expected to take much longer than inside the TS, probably even longer 
than 11 or 22 years of the solar activity- or magnetic-cycles. In this case, the plasma regions 
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originating from the active and quiet sun are both expected to exist alte rnatively in the 



20071 ) (Dr. H. Washimi, 



heliotail, as confirmed by a recent MHD simulation (I Washimi et al. 
2009, private communication). If the observed attenuation of the sidereal anisotropy occurs 
predominantly in the heliotail, therefore, it is possible to expect the attenuation magnitude 
showing no strong and direct correlation with the solar activity- and/or magnetic-cycles, as 
we saw in the observation with Matsushiro. This is just one of possible interpretations and 
we need to confirm it by analyzing the GCR propagation in the model heliosphere including 
the heliotail predicted by MHD simulations. 

Finally we note that the magnetic structure of the heliosheath is largely unknown. The 
possible role of the heliosheath is particularly interesting in the light of the recent finding of 
the IBE X spacecraft, mapping t he global heliosphere by detecting energetic neutral atoms 



(ENAs) flMcComas et al. 



2009i ). The first results of IBEX showed that ENAs are sensitive 



to, and give information on the struct ure of the large-scale magnetic field surrounding the 



heliosphere (IHeerikhuisen et al. 



20101 ). TeV cosmic rays also sense the global structure 
of the magnetic field within the whole heliosphere as well as the field surrounding the 
heliosphere, and can provide an additional tool to explore the global structure of these 
fields. 



The continuous observation with the Matsushiro underground muon detector has been 
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Table 1. Amplitude and phase of the yearly mean diurnal anisotropy observed by 

Matsushiro in 1985-2008. 



si 



so 



ASI 



year 



a? (%) 



4* (h) 



of (%) 



4$ (h) 



of (%) 



»f (h) 



1985 


0.045±0.013 


2.66±2.00 


0.022±0.009 


8.80±1.58 


0.009±0.009 


23.53±4.04 


1986 


0.034±0.012 


4.57±1.82 


0.024±0.009 


6.16±1.40 


0.011±0.009 


23.09±2.96 


1987 


0.032±0.013 


23.53±0.50 


0.035±0.009 


6.70±1.01 


0.003±0.009 


7.10±13.58 


1988 


0.069±0.012 


1.38±1.29 


0.042±0.009 


5.84±0.80 


0.013±0.009 


16.34±2.56 


1989 


0.022±0.012 


23.66±0.39 


0.009±0.009 


9.42±3.89 


0.019±0.009 


6.02±1.82 


1990 


0.051±0.013 


0.95±0.94 


0.024±0.009 


8.81±1.50 


0.012±0.009 


10.13±2.96 


1991 


0.065±0.013 


3.98±1.79 


0.016±0.009 


11.70±2.29 


0.025±0.009 


22.49±1.44 


1992 


0.055±0.012 


2.96±2.00 


0.029±0.009 


11.36±1.21 


0.014±0.009 


15.04±2.45 


1993 


0.047±0.013 


6.01±1.03 


0.008±0.009 


6.57±4.54 


0.022±0.009 


2.16±1.62 


1994 


0.030±0.013 


21.77±1.97 


0.043±0.009 


10.41±0.84 


0.028±0.009 


7.96±1.28 


1995 


0.047±0.013 


0.13±0.13 


0.028±0.009 


8.09±1.25 


0.002±0.009 


12.97±15.95 


1996 


0.010±0.013 


21.46±3.58 


0.038±0.009 


8.70±0.93 


0.013±0.009 


8.27±2.75 


1997 


0.028±0.013 


1.95±1.83 


0.046±0.009 


5.22±0.75 


0.012±0.009 


7.09±2.96 


1998 


0.013±0.013 


22.64±1.82 


0.040±0.009 


5.92±0.87 


0.018±0.009 


13.05±1.93 


1999 


0.029±0.013 


3.09±2.25 


0.015±0.009 


4.23±2.34 


0.022±0.009 


2.68±1.58 


2000 


0.048±0.013 


2.06±1.76 


0.021±0.009 


8.01±1.70 


0.016±0.009 


18.33±2.26 


2001 


0.052±0.013 


4.51±1.60 


0.033±0.009 


8.34±1.06 


0.030±0.009 


0.05±1.19 


2002 


0.053±0.013 


2.15±1.79 


0.034±0.009 


10.76±1.03 


0.021±0.009 


20.12±1.64 


2003 


0.046±0.013 


4.35±1.75 


0.027±0.009 


12.95±1.30 


0.033±0.009 


22.32±1.08 


2004 


0.023±0.013 


4.28±2.52 


0.023±0.010 


11.81±1.61 


0.007±0.010 


3.22±5.23 


2005 


0.020±0.014 


23.98±0.03 


0.050±0.010 


7.78±0.79 


0.013±0.010 


7.21±3.09 


2006 


0.035±0.013 


4.39±1.91 


0.009±0.009 


6.41±3.79 


0.008±0.009 


22.15±4.20 


2007 


0.034±0.013 


2.48±2.03 


0.033±0.009 


6.07±1.07 


0.009±0.009 


23.49±3.85 


2008 


0.033±0.013 


2.48±2.03 


0.028±0.009 


8.84±1.30 


0.003±0.009 


1.76±10.31 



a Errors arc deduced from the muon counts. The amplitude and phase </>^ in SI time frame 
are corrected for the "sideband" effect (see text). 
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Fig. 1. — Summation dials of the yearly mean harmonic vector observed by Matsushiro 
in 1985-2008. Yearly mean diurnal vectors are displayed in the harmonic dial, where the 
amplitude (af ) of each yearly vector is represented by its length from the origin, while the 
phase (0f ) is represented by the angle measured clockwise from the vertical (+y) axis, i.e. 
the local times of +y, +x, —y, and —x axis directions in this figure are 00:00, 06:00, 12:00 
and 18:00 in the local time (LT), respectively. To demonstrate the long-term variation, we 
plot "vector summation dials" in this figure where yearly vectors are summed one by one 
from the first year of 1985. We also divide each yearly vector by 24, the number of years of 
observation. The final data point in each summation dial, therefore, represents the average 
yearly vector. The average vectors in the SI, SO and ASI time frames are shown by full circle, 
full square and open triangle, respectively. The vectors in the SI time frame are corrected 
for the "sideband" effect (see text). The uncorrected average SI vector is also displayed by 
a gray circle. Error bars of each average vector are deduced from the dispersion of x- and 
^/-components of 24 yearly mean vectors. 





1995 Year 2000 



Fig. 2. — Long-term variation of the diurnal anisotropy observed by Matsushiro in 1985- 
2008. Yearly mean amplitude (solid circles) and phase (open circles) of the solar diurnal 
anisotropy (b) and the sidereal diurnal anisotropy (c) are plotted as functions of year. The 
sidereal diurnal anisotropy in (c) is corrected for the "sideband" effect (see text). The yearly 
averages of the sunspot number (SSN) and the neutral sheet tilt angle (TA) representing the 
solar activity are also shown in panel (a). Plotted in panel (b) are the amplitude and phase 
of the "additional" anisotropy which is derived by subtracting the harmonic vector expected 
from the Compton- Getting anisotropy (with an amplitude of 0.035 % and a phase of 05:56 
LT) from the observed anisotropy. The yearly mean phases in years when the amplitude 
is below 0.02 % are omitted from the plot in order to exclude large fluctuation due to the 
insignificant anisotropy. Errors in panels (b) and (c) are statistical errors deduced from the 
muon counts. The gray solid curves in (b) and (c) also show the 12-month central moving 
average of the amplitude. The yearly mean SSN and TA are calculated respectively from the 
monthly SSN value ( http: / /solarscience.msfc.nasa.gov[ ) and the Carrington rotation average 
of the TA computed using a model based on the radial boundary condition at the photosphere 
(http://wso.stanford.edu/Tilts.html). Errors in panel (a) are deduced from dispersions of 
the monthly mean SSN and the Carrington rotation average of TA used for calculating 
yearly mean values. The periods of the solar polar magnetic field reversals deduced from the 
solar polar magnetic field data (http://wso.stanford.edu/Polar.html) are also indicated by 
the vertical shaded epochs. 
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Fig. 3. — Average solar daily variations R*(ti) in two periods selected according to the 
amplitude of the additional anisotropy. Panel (a) shows by solid circles the average R*(ti) in 
6 years of 1991, 1992, 1994 and 2002-2004 (period I) when the amplitude of the "additional" 
anisotropy in Fig. 2b exceeds 0.035 %, while panel (b) displays R*(U) in 6 years of 1986- 
1988, 1997, 1998 and 2007 (period II) when the amplitude is below 0.015 %. Dashed curve in 
panel (b) represents the diurnal variation expected from the CG anisotropy arising from the 
earth's orbital motion around the sun. In panels (a) and (b), R*(ti)s are shown as functions 
of the local solar time ti on the horizontal axis. Error is deduced from the dispersion of the 
yearly mean of R*(ti) used for the calculation of the average and displayed by gray curves. 



(a) period I 




(b) period II 
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Fig. 4. — Average sidereal daily variations R*(ti) in four periods selected according to the 
solar activity- and magnetic-cycles. All daily variations in this figure are corrected for the 
"sideband" effect (see text). Panels (a) and (b) display by solid circles R*(U)s in "active" 
(1988-1992, 1999-2002) and "quiet" (1985-1987, 1994-1997, 2005-2008) periods of the solar 
activity, respectively, while panels (c) and (d) show respectively R*(U)s in "positive" (1991- 
1998) and "negative" (1985-1988, 2001-2008) periods of the solar polar magnetic field. Error 
is deduced from the dispersion of the yearly mean of R*(ti) used for the calculation of each 
average and displayed by gray curves. 
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Fig. 5. — The Single Band Valley Depths (SBVDs) observed by Matsushiro and Milagro 
experiment in a 7-year period between 2000-2007. The SBVDs by Matsushiro and Milagro 
experiment are plotted by full and open circles, respectively, as functions of year on the 
horizontal axis. The SBVD by Matsushiro in this figure is corrected for the "sideband" 
effect and multiplied by three in order to compensate for the attenuation of the amplitude 
in sub-TeV region. The yearly mean SBVD by Matsushiro is calculated from 12-month data 
between June and July adjusting to the analysis period by Milagro experiment (Abdo et 
al., 2009). The gray solid curve also displays the 12-month central moving average of the 
SBVD by Matsushiro. The error of Matsushiro's SBVD which is corrected for the "sideband" 
effect is only statistical, while the error of Milagro's SBVD includes both the statistical and 
systematic errors evaluated from the amplitude of the observed ASI diurnal variation (see 



Abdo et al. 



20091 ). 



